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Abstract 

The CACTUS atmospheric Cherenkov telescope collaboration recently reported a gamma-ray excess 
from the Draco dwarf spheroidal galaxy. Draco features a very low gas content and a large mass-to- 
light ratio, suggesting as a possible explanation annihilation of weakly interacting massive particles 
(WIMPs) in the Draco dark-matter halo. We show that with improved angular resolution, future 
measurements can determine whether the halo is cored or cuspy, as well as its scale radius. We find 
the relevant WIMP masses and annihilation cross sections and show that supersymmetric models can 
account for the required gamma-ray flux. The annihilation cross section range is found to be not com- 
patible with a standard thermal relic dark-matter production. We compute for these supersymmetric 
models the resulting Draco gamma-ray flux in the GLAST energy range and the rates for direct neu- 
tralino detection and for the flux of neutrinos from neutralino annihilation in the Sun. We also discuss 
the possibility that the bulk of the signal detected by CACTUS comes from direct WIMP annihilation 
to two photons and point out that a decaying-dark-matter scenario for Draco is not compatible with 
the gamma-ray flux from the Galactic center and in the diffuse gamma-ray background. 
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1 Introduction 

Despite compelling indirect astrophysical and cosmological evidence, the fundamental nature of non- 
baryonic dark matter remains elusive (see Refs. P3E] f° r recent reviews). Upgrades of direct detectors 
looking for scattering of dark-matter particles from nuclei, future neutrino telescopes, and space- 
based antimatter and gamma-ray detectors will dramatically enhance the possibility to discover dark- 
matter particles in the Galactic halo. Given our complete ignorance on the nature of dark matter, 
any anomalous experimental result that might point to dark matter should be carefully analyzed, 
and possibly cross-correlated both with other observational constraints and with existing theoretical 
models. 

Among indirect-detection techniques that seek products of dark-matter annihilations, those search- 
ing for neutrinos and photons play a special role. Unlike electrically charged particles, like antiprotons 
or positrons, which diffuse in the Galactic magnetic field, the arrival direction of neutrinos and pho- 
tons points to where the dark-matter annihilation took place. In particular, the forthcoming launch 
of the Gamma Ray Large Area Space Telescope (GLAST) [3] and the rapidly developing field of 
ground based Atmospheric Cherenkov Telescopes (ACTs) 0] make the search for energetic photons 
from dark- matter annihilation a particularly promising and exciting endeavor. 
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Since the rate per unit volume for dark-matter annihilation depends on the square of the dark- 
matter density, there may be great advantage to seeking astrophysical locations where the dark-matter 
density is believed to be high. The center of the Milky Way has been viewed as a promising target as 
it might host a dark-matter spike and/or because the Galactic halo could feature a steep cusp towards 
its center jHj. On the other hand, diffuse gamma-ray backgrounds produced by the spallation of 
cosmic rays on interstellar gas, and the uncertainties on the distribution of molecular hydrogen in the 
galactic ridge, may blur the dark-matter- induced signal (but for a different viewpoint see Ref. 0). The 
presence of other gamma-ray emitters in the central region of our Galaxy, including those associated 
with the supermassive central black hole and the supernova remnant Sgr A*, however, also make it 
difficult to distinguish a putative dark- matter signal and its astrophysical background [Tj |S| . 

It has therefore been suggested that one should instead seek dark-matter-induced gamma rays 
from other sources that feature large dark-matter densities but that remain devoid of background 
sources. In particular, several extragalactic sources have been considered (§1 E3 EEU EE21 E3- Among 
these, dwarf spheroidal (dSph) galaxies are among the most dark-matter dominated systems, featuring 
mass-to-light ratios as high as M/L ~ 250Mq/Lq |14j . Moreover, at least four dSph galaxies with 
very large M/L lie within 100 kpc from the Milky Way center. 

The dark-matter-induced photon flux depends on the dark-matter distribution, and particularly on 
its inner structure, which is unfortunately poorly known in the case of local dSphs. Ref. ^U] modeled 
the dark-matter halos of the four nearby dSph, Carina, Draco, Ursa Minor, and Sextans, with King 
profiles. Even with optimistic assumptions on the background-rejection capabilities of ACTs, they 
found that the gamma-ray flux expected from supersymmetric models from local dSphs is typically at 
least three orders of magnitude below the background. They find that if the dark-matter distribution 
is clumpy, the signal can be boosted by at most a factor of 40, which would still be insufficient. A 
re-examination of the Draco dSph galaxy in Ref. jHj arrived at somewhat different conclusions. This 
work showed that, assuming a steeply cusped isothermal dark-matter halo, a large portion of the 
supersymmetric parameter space could produce a signal visible at both forthcoming ACTs and at 
GLAST. This work also pointed out that, depending on the magnetic-field strength inside Draco, the 
radio-continuum limits for Draco obtained in Ref. |15j could also rule out sizeable portions of the same 
parameter space. Ref. jllj conducted an extensive analysis of the 44 nearest galaxies in the Local 
Group, pointing out other promising extragalactic candidates, such as M31. A systematic comparison 
between the case of the Galactic center and that of various nearby dSph galaxies (Sagittarius, Draco, 
and Canis Major) was carried out in Ref. ^3], where a wide array of halo profiles was also employed. 
Depending on the latter, Ref. ^3] concludes that if dark matter is supersymmetric, dSph galaxies may 
indeed produce a detectable signal in ACTs and at GLAST. 

Recently, a gamma-ray excess from the direction of Draco has been detected by the Solar 2 Heliostat 
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Array CACTUS, located in Barstow, California Although the robustness of the signal has still 
to be tested, the possibility of ascribing the excess over the off-source background to dark-matter 
annihilation in Draco's halo is certainly intriguing, as Draco is a dark-matter dominated system (see 
Sec.[2J) that is not expected to host any other significant gamma-ray source This possibility was 

recently envisaged in Ref. ^H]. In the present analysis, we study the impact on the gamma-ray flux 
of various, astrophysically motivated, halo models ^HI- We point out that if the angular resolution 
of the experiment can be improved and understood, then measurement of the angular distribution of 
the gamma-ray excess can discriminate between cuspy (e.g., the Navarro- Frenk- White profile 20 ) and 
cored profiles and determine the halo scale radius (Sec. 13.1]) . In the absence of a full analysis of the raw 
counts reported by CACTUS [T2|, the estimate of the flux of photons to be attributed to dark-matter 
annihilation appears particularly critical. Here, we conservatively estimate the putative gamma-ray 
flux detected by CACTUS ^Hj bracketing it between an upper estimate where all the CACTUS 
excess counts over the background are attributed to the dark-matter signal (following the approach of 
Ref. 18 J, and a lower estimate where only the counts in the innermost angular region are supposed 
to originate from dark-matter annihilations. We then scrutinize the dark-matter interpretation of the 
CACTUS excess from Draco on model-independent grounds on the planes defined by the particle mass 
and annihilation cross section (Sec. 13.21) . We specialize to the particular case of supersymmetric dark 
matter in Sec. 13.31 where we also make predictions for the expected detection rates for models consistent 
with the Draco excess with other dark-matter searches. We consider in Sec. l3.4l the possibility that the 
bulk of the excess photons detected by CACTUS come from the monochromatic 77 line. In Sec. 13.51 
we show that the signal cannot come from the decay of a very long-lived dark species. We finally draw 
our conclusions in Sec. 0] 

2 The dark-matter halo of Draco 

Draco was the first dSph galaxy to show evidence of a large dark-matter content. A total mass 
much larger than the amount of visible matter was inferred in Ref. [2j| from the measurement of four 
carbon stars in Draco as early as 1983. A few years later, the assessment of the velocities of 15 more 
stars allowed Ref. |22| to compute a mass-to- light ratio larger than 50 Mq/Lq, a result subsequently 
confirmed by a much larger sample of 91 Draco stars in Ref. [2S] and 17 stars in Ref. [21]. More 
recently, the existence of an extended dark-matter halo was inferred in Refs. |25[ I26j from a sample 
of radial velocities of 159 giant stars out to large projected radii; the previously assumed hypothesis 
that mass follows light was shown to be inconsistent at the 2.5-<t confidence level, and tentative 
evidence for a nearly isothermal dark-matter distribution was provided. A new re-evaluation of the 
Draco dark-matter distribution, based on new data [27j and on two-component high-resolution N- 
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body simulations, together with cosmological predictions for the properties of dark-matter halos, lead 
Ref. 19 to a few important conclusions, relevant for the present discussion. First, it was pointed 
out that both a cored and a cuspy profile are compatible with the observational data on Draco and 
with the results of numerical simulations. Second, the possibility that Draco is a tidal dwarf (i.e., a 
virtually unbound stellar stream tidally disrupted in the Milky Way gravitational potential |28| ) was 
ruled out. This might very well be in contrast with other local dwarf galaxies |19l I13j . for which a 
smaller amount of data is available. 

The analysis of Ref. [TO] assumes two types of (spherical) dark-matter halos. The first one, moti- 
vated by iV-body simulations, is the well known Navarro- Frenk- White (NFW) density profile |20| . 



PNFw( r ) — 7 / rrz ; 7 (1) 
(r/r s )(l + r/r s y 

as a function of radius r, while the second is the observationally motivated Burkert profile [2 



PBur(r) (l + r/r s )[l + (r/r s yy (2) 

The flux of gamma rays induced by dark-matter annihilation is clearly very sensitive to the parameters 
entering the two profiles; i.e. the scaling density (p s ) and scale radius (r s ). The (p s ,r s ) plane is 
constrained by observations as well as cosmological arguments. We reproduce in Fig. ^ with gray 
shading the viable portions on the (p s , r s ) plane determined in Ref. ^2] • The leftmost region is 
ruled out by the requirement that the dark-matter halo formed early enough to allow the subsequent 
formation of the bulk of its stellar component. The upper-right limit results from the requirement that 
the number of dSph galaxies in the Local group is equal at least to the number of observed objects. 
Finally, the viable region is bounded from below by the condition that Draco's virial radius extends 
out at least to the last observed point in its surface-brightness profile and by consistency with the 
observed line-of-sight stellar velocity dispersion. We indicate with blue crosses points corresponding to 
models that have \ 2 < 9.5 between the modeled and observed line-of-sight velocity-dispersion profiles 
according to the analysis of Ref. |19j . For future convenience, we pick four benchmark halo models, 
among those analyzed in Ref. that feature a wide range of p s and r s ; we highlight these models 
in Fig. ^ with red squares. We collect the scaling densities and radii, as well as the values of J_3 and 
J_5 for the four benchmark models, in Tabled 

The observed integral flux of gamma rays with energies E > E~ from dark-matter annihilation per 
unit solid angle is 

iE~, _ NJav) 
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Jaq, (3) 

'DM 

and depends on the particle mass todM; on its annihilation cross section (crv), on the number iV 7 of 
photons with energies E > E~ produced per annihilation, and on the integral Jac along the line of 
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Figure 1: Curves of constant line- of- sight integral of the square of the dark-matter density Jasi (see Eq. 
for the definition of this quantity), for Aft — 10~ 5 (solid black lines) and Afl — 10 -3 (dashed red lines), on 
the (p s ,r s ) plane. For an NFW (left) and for a Burkert (right) dark-matter halo, gray regions correspond to 
cosmologically allowed ranges of the two halo parameters, while models indicated with blue crosses also give 
consistent simulated stellar velocity dispersions. We highlight with red squares the four benchmark halo models 
of Tabled 



sight of the square of the dark- matter density, averaged over the solid angle AO of the detector, 
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(5) 
(6) 



so = 75.8 ± 0.7 ± 5.4 kpc being Draco's heliocentric distance, and rt its tidal radius. The latter was 
evaluated, following the approach of Ref. JH]; according to the Roche criterion (the radius such that 
the average mass in the dSph is equal to the average interior mass in the Milky Way halo); different 
tidal radii correspond to different Milky Way (and, though more weakly, to different Draco) dark- 
matter halos, ranging within less than one order of magnitude. However, since Jac depends quite 
weakly on rt^>r s , we resorted to an isothermal profile for the Milky Way, which typically gives rt ~ 7 
kpc. 
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Table 1: Input parameters and line-of-sight integrals for four benchmark Draco halo models from Ref. 

We show in Fig. ^ iso-level curves for the quantities J_3 = </a£2=io- 3 sr an d J-5, relevant, respec- 
tively, for the full one-degree angular region where the gamma-ray excess from Draco was observed 
and for the innermost 0.1-degree angular region, where the largest dark- matter-induced gamma-ray 
flux is expected, on the (p s ,r s ) plane, for the NFW profile (left) and for the Burkert profile (right). 
We define here and the remainder of the paper the quantities Jac in units of 10 23 GeV 2 cm -5 . We 
note that (1) the observationally and cosmologically consistent range for J_5 spans roughly one order 
of magnitude in the case of a cored profile and 2 orders of magnitude for a cuspy profile. (2) The 
ratio J-3/J-5 ~ 1 grows significantly with r s for cored profiles, while J-3/J-5 ~ 10 4- 100 for cuspy 
profiles, with smaller values corresponding to larger scale radii. We find as a consistent range for J_5 
(again in units of 10 23 GeV 2 cm" 5 ), 

0.11 < J_5 < 2.96, Burkert profile, 
0.35 < J_ 5 < 64.3, NFW profile. 

3 The CACTUS excess and dark matter 

The CACTUS ACT observed a gamma-ray excess over background from an angular region extending 
approximately 1 degree around the direction of the Draco dwarf spheroidal galaxy ^B]. The CACTUS 
energy threshold is around 50 GeV, and no statistically significant excess for energies greater than 150 
GeV was reported. 

Given several potential issues concerning the effect of the integrated starlight from the back- 
ground and from Draco's stars, the noise-reduction procedures, and the intrinsic background related 
to misidentified electromagnetic showers induced by hadrons and electrons, it is extremely difficult 
to reliably estimate the signal gamma-ray flux from Draco. Some portion of the excess counts over 
the background reported by CACTUS will presumably not be related to the putative dark-matter 
signal. It seems reasonable to assume that the gamma-ray flux is bracketed between two extremes: as 
a conservative upper limit, one can consider the overall excess counts detected from the angular region 
of 1 degree centered around Draco above the average background measured outside Draco, which cor- 
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responds to approximately 30,000 photons, detected with an effective area of the order of 5 x 10 4 m 2 
in 7 hours of observation Hi . (This procedure was applied in the recent analysis of Ref. |18j.) 

As an alternative, we can proceed as if the experiment had an angular resolution of 0.1 degrees. 
Currently, the angular resolution toward Draco is not yet well understood by the CACTUS collabora- 
tion. A point-source resolution of 0.3 degrees was obtained toward the Crab nebula, but the resolution 
toward Draco may be even poorer. It is thus premature to attribute the 1 degree spread in the Draco 
excess to a source of a 1-degree spatial extent; it may well still be consistent with a point source. Still, 
to illustrate the possibilities with future measurements with better resolution, we proceed with our 
theoretical investigation as if the experiment had such a resolution. In this illustrative exercise, one 
can then suppose that the signal comes only from the innermost 0.1 degree angular region around the 
center of Draco, and that the rest of the excess is due to spurious effects. This second procedure is 
found to be equivalent, in the estimate of the signal flux, to requiring that in the innermost region of 
Draco, which should contain most of the dark-matter-induced gamma-ray flux, the ratio of the signal 
over the estimated off-source background for an ACT like CACTUS reproduces the observed relative 
excess size 20%) [TB] . 

An angular region AQ = 10 -3 sr, corresponding to an angular radius of approximately 1 degree, 
gives a signal flux (above the CACTUS energy threshold of 50 GeV) of <f)% = <f> 1 (E 1 > 50 GeV) sa 
2.4 ~r~ 3.4 x 10~ 9 cm~ 2 s~ 1 , where the lower value corresponds to an effective area A e g = 5 x 10 4 m 2 , 
while the upper value to an average effective area which takes into account the energy dependence 
over the interval 50 < £ 7 /GeV < 250 pH EE] • 

In the second approach, based on the innermost 0. 1-degree angular region to which an ACT may 
ultimately be sensitive (corresponding to AQ ~ 10 -5 sr), the main sources of background correspond 
to misidentified gamma-like hadronic showers and cosmic-ray electrons. The diffuse gamma-ray back- 
ground should instead not contribute significantly, but we will take it into account as well. We use 
the following estimates for the ACT background [3j: 



where eh a d 1 parameterizes the efficiency of hadronic rejection, and where we took for the diffuse 
gamma-ray background the most conservative differential spectral index (corresponding to the extra- 
galactic gamma-ray background) and normalized the flux to the value of the diffuse emission above 1 
GeV from EGRET [SO] in the direction of Draco (/ = 86.4°, b = 34.7°). Although the Galactic diffuse 






diVdiff 




(9) 
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background typically exceeds the extragalactic background at EGRET energies, the two spectra are 
such that the extragalactic background should dominate above 50 GeV. We therefore derive an overall 
background level, using an angular acceptance of 10 _5 sr appropriate for ACTs under our assumptions, 
of 

& g = 0bck g (£ > 50 GeV) ~ 1.2 x 10~ 10 cm" 2 s~\ (10) 
<j>b£ g = 4>bck g (E > 150 GeV) ~ 1.5 x 10~ n cm^s" 1 . (11) 

Comparing the signal outside Draco and towards its center, the observed excess is around 20% of 
the background, which gives a putative gamma-ray excess 0^° 5 « 2.4 x 10 -11 cm~ 2 s _1 . We let here the 
signal vary within one order of magnitude around that central value, to account for the uncertainties 
in the background estimation and in the actual size of the claimed excess. We therefore consider a 
signal range, in this approach, 0.06 < ~ 0.6. As mentioned above, this second conservative 

estimate of the signal is consistent with the actual number of counts reported by CACTUS within an 
angular radius of around 0.1 degrees. 

No excess flux has been observed from Draco above 150 GeV. This leads to a further constraint 
on the dark-matter interpretation. We consider two putative upper limits: the strongest comes from 
the requirement that in the central bins the signal is less than 5% of the ACT backgrounds, and reads 
0-5° ^ 7.5 x 10~ 13 cm _2 s _1 ; the most conservative requirement is instead that the signal flux does 
not exceed the Poisson fluctuation of the actual number of counts above the energy threshold of 150 
GeV; this gives a limit on the signal flux of (j) 15 ® < 12.6 x 10 -12 cm _2 s _1 . 

In the following, we will refer to the CACTUS signal as a photon flux in the range 

4>% < 3.4 x 10~ 9 cm-V 1 (f>% > 7.2 x 10~ 12 cm^V 1 <^ 5 5 ° < 0.75 -=- 12.6 x 10~ 12 cm -2 s -1 . (12) 

Again, the AQ = 10~ 5 numbers do not describe the current Draco data; rather, they describe results 
of a hypothetical experiment with 0.1 degree resolution that look like the current CACTUS results. 

3.1 The gamma-ray angular distribution 

The gamma-ray excess is spread over 1 degree. However, given uncertainties in the current CACTUS 
angular resolution around Draco, we cannot currently use the observed spread to discriminate between 
a source with a 1 degree spatial extent and a point source. Still, it is conceivable that forthcoming 
ACT measurements may achieve a resolution as good as 0.1 degrees. In this Section, we illustrate how 
future such measurements may shed light on the halo profile for Draco. To do so, we proceed as if the 
resolution of CACTUS were in fact already 0.1 degrees. 

Our second estimate of the flux, depending on the assumed dark-matter distribution, and particu- 
larly on the scale radius of the dark-matter halo, typically gives a total flux from the central 1 degree 
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Figure 2: A comparison among theoretical predictions for the angular distribution of the gamma-ray flux from 
dark-matter annihilation and the excess reported by CACTUS for four different Draco dark-matter halos \iyjj . 
featuring a wide range of scale radii, and giving consistent N -body-simulated stellar-velocity dispersions. Keep 
in mind that the current CACTUS data may have a resolution of roughly one degree, and so conclusions about 
the Draco halo cannot yet be drawn by comparing the current data with the theoretical curves. 

angular region which can be significantly lower than the total counts reported by CACTUS. Within 
this approach, the counts originating from regions outside the center of Draco would correspond to 
photons produced by dark-matter annihilation only if a very large scale radius (of the order of 1 kpc) 
were assumed. 

What we will now show is that if CACTUS excess is real and if the Draco halo is cuspy, then an 
ACT measurement with an 0.1 degree resolution should see a very peaked angular distribution. To 
do so, we reproduce in Fig. [2] the measured photon excess in the 1-degree angular region centered 
on the location of Draco. We compare the putative excess with the predictions, for the photon-flux 
angular distribution, stemming from four different dark-matter halos. We normalize the photon fluxes 
for the various halos to the total number of excess gamma rays reported by the CACTUS preliminary 
results. The four profiles were chosen among those quoted in Table 1 of Ref. |19j : two of them are 
Burkert profiles with scale radii of 1.4 kpc (BUR1, profile Bl in Ref. and 0.35 kpc (BUR2, profile 
B2 in Ref. ^H])) and two of them are NFW profiles with scale radii of 2.8 kpc (NFW1, profile Nl in 
Ref. [Ej) and 0.18 kpc (NFW2, profile N5 in Ref. [El). We deduce from Fig. that if the current 
resolution of CACTUS were as good as 0.1 degrees, then the NFW profiles would produce an excess 
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of gamma rays in the central bins over what is now seen. If, however, the current angular distribution 
is still observed even with an 0.1 degree resolution, then cored profiles, with scale radii of order 1 kpc, 
would be indicated by the data. We conclude therefore that improved angular resolution is warranted 
to understand better the halo structure. 

3.2 Implications for particle dark matter 

No gamma-ray source in the direction of Draco was identified by EGRET in its all-sky survey. An 
analysis was carried out by the EGRET collaboration in Ref. jHJ. The absence of a point-like source 
from the direction of Draco can be independently used to draw an upper limit on the integral flux of 
photons above a threshold of 1 GeV, following the analysis of Ref. |S2]- As in Ref. [5], one can take, 
as the EGRET upper limit on point-like sources gamma-ray fluxes the flux of the least significant 
point-source detection in the EGRET catalog, corresponding to the Large Magellanic Cloud, which 
translates into the requirement 

0l 3 = </>(E>l GeV) < Regret - 10~ 8 cm-V 1 . (13) 

This flux agrees with the theoretical estimate for the flux sensitivity of EGRET to point sources 
determined in Ref. [3*3*] . 

The actual EGRET data from the direction of Draco were collected, in Ref. [HJ, in 7 energy bins, 
featuring different angular cuts and different exposures. Every energy bin is also accompanied by 
a background estimate. In order to compare the aforementioned point-source sensitivity of EGRET 
with the actual photon count, one needs to pick a halo profile for Draco and a signal photon spectrum. 
We also clustered the four lowest-energy bins (0.1 < E^/GeV < 1) and the three highest-energy bins 
(1 < -E-y/GeV < 10), since in some bins the background estimate exceeds the photon count, and in the 
highest-energy bins the statistics are very poor. We require that in the two energy intervals the signal 
does not exceed the difference between the measured photon counts and the estimated background. 

We find that in general the point-source sensitivity agrees, within a factor at most 2, with the 
bound from the highest-energy bins, while that from the smallest-energy bins is approximately one 
order of magnitude weaker (and, in particular, the harder the signal photon spectrum, the weaker the 
latter bound). The bound from the highest-energy bins tends to be stronger than the point-source- 
sensitivity constraint for (1) smaller particle masses, (2) for cored rather than cuspy dark-matter 
halos, and (3) for softer photon spectra. Since, however, the point-source-sensitivity criterion appears 
to be less dependent on the poor statistics of the highest-energy bins (consisting of only a few photon 
counts), and since the two criteria essentially agree, we hereafter indicate as the EGRET bound the 
limit in Eq. (fTSj). 
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WIMP mass [GeV] WIMP mass [GeV] 



Figure 3: Regions on the (m WIMP , (crv)) plane compatible with the CACTUS excess and with other observations 
for a soft photon spectrum (of the same type as that shown with red lines in Fig. [3| bb final state). The two left 
panels assume the cuspy NFW2 dark-matter halo for Draco, while those on the right the cored BUR1 halo (see 
Tabled . 

The VERITAS collaboration also observed Draco with the Whipple-lOm ACT, reporting a null 
search for high-energy gamma rays above a threshold E~ ~ 400 GeV, which translates into the bound 

<f> m ! = <f>(E > 400 GeV) < <^ ipple _ 10m * 1 x lO" 12 corV 1 . (14) 

The number of gamma rays from the annihilation of a WIMP integrated over cin energy E^ Ccin 
be written as 



N^E,) = 




where the symbol / refers to any WIMP-annihilation final state, yielding a gamma-ray spectral func- 
tion (differential number of photons per WIMP annihilation) dNy/dE. Different final states give 
different spectral functions. We take here as benchmark cases those giving the hardest and the softest 
spectra among the final states that are relevant in the case of supersymmetric dark matter; i.e., t + t~ 
and bb, respectively. 

We show in Figs. El and 0] the region, on the (m wlMP , (av)) plane, compatible with the putative 
CACTUS signal quoted in Eq. I|12|). and the various constraints discussed above, including the two 
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Figure 4: As in Fig. but for a hard photon spectrum (blue dashed lines in Fig. [3| t + t final state). 

estimates of the limit from the non-observation of a significant excess of gamma rays with energies 
above 150 GeV. Fig. prefers to a soft photon spectrum, corresponding to the final state / = bb, while 
Fig. 01 to a hard photon spectrum, / = t + t~. The two panels to the left to a cuspy profile (namely, 
the NFW2 profile, i.e. profile N5 in Ref. JH]), while the two panels to the right employ a cored 
dark-matter halo (namely, the BUR1 profile; i.e., profile B3 in Ref. |19|). 

We emphasise how the range of dark-matter annihilation cross sections needed to explain the 
CACTUS signal is largely above what a naive estimate from the thermal production of relic WIMPs 
would suggest. Assuming that the pair annihilation rate (cry) is energy independent, and there- 
fore that annihilations proceed through s-wave processes, the WIMP \ relic abundance scales as 
£l x h? ~ 3 x 10~ 27 cm s / (crv) Requiring that all the CACTUS excess originates from dark-matter 
annihilations would therefore naively entail thermal relic abundances as low as ^ x /i 2 ~ 10 -5 -j- 10~ 6 . 
For reference, we indicate the {av) range deduced from the previous relation and from the 2-cx WMAP 
estimate of the CDM abundance with light blue bands in fig. and |IJ Resonant annihilations, 
coannihilation processes or thresholds can significantly distort, though, the above outlined estimate: 
we postpone to sec. 13.31 (see in particular fig. |HJ) a discussion of the compatibility of the CACTUS 
signal with the thermal production of dark matter for the particular case of neutralinos. 

We also include the sensitivity of GLAST, computed assuming a total exposure of 3.2 x 10 11 cm 2 s 
and the diffuse gamma-ray background measured by EGRET, Eq. Q, and requiring the strongest 
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among the following conditions: (1) Ns/y/Ns > 5, N$ and Nb being the total number of signal and 
background events, or (2) N$ > 10. The angular acceptance was set to maximize the signal-to-noise 
ratio (namely, the quantity vAO for a given halo profile) between the minimal angular resolution 
and the maximal field of view. 

In the case of a soft photon spectrum (Fig. the EGRET bound rules out large parts of the 
CACTUS signal, regardless of the assumed halo profile. With a soft photon spectrum, all models 
consistent with the lower estimate of the CACTUS signal will be largely within the GLAST sensitivity. 

With a hard signal photon spectrum (Fig. |1J), the bound from the Whipple-lOm observation of 
Draco is significantly more effective than with a soft spectrum; again, almost all the region compatible 
with the CACTUS signal will be within GLAST reach. We also remark that with the BUR1 profile 
(Fig. right), and with a hard photon spectrum, we find models that (1) give a total photon flux 
equal to the CACTUS counts; (2) give an angular distribution of the photon counts that is very close 
to that reported by CACTUS (see Fig. 12 keep in mind that the theoretical curves reported in Fig[2] 
refer to a much better angular resolution than that of CACTUS); and (3) are consistent with the 
conservative criterion on the gamma-ray flux above 150 GeV. 

The two black crosses indicate two benchmark points, whose differential photon spectrum is re- 
produced in Fig. El for illustrative purposes, together with the expected diffuse gamma-ray back- 
ground, for a space-based gamma-ray telescope with angular resolution AS7 = 10~ 3 sr (EGRET-like) 
and Af2 = 10 -5 sr (GL AST-like). We also show the monochromatic gamma-ray line, assuming 
{av) 17 /(av)t t = 3 x 10~ 4 (close to the maximal value for super symmetric dark matter; see Sec. 13.411 . 
and different energy resolutions, appropriate for the two space-based detectors. 

3.3 Supersymmetric dark matter 

Our analysis was based, up to this point, on a model-independent approach as far as the annihilation 
cross section and the final-state branching ratio for dark-matter annihilation are concerned. In order 
to predict the rates for dark-matter detection in other search arenas, such as direct detectors and km 3 
neutrino telescopes, one needs to consider specific particle models. To this extent, we now specialize 
to supersymmetric dark matter and consider the minimal CP-conserving supersymmetric extension of 
the standard model (MSSM), and perform a random scan over its parameter space. For all models, we 
impose constraints from direct supersymmetric-particles searches at accelerators, rare processes with a 
sizeable potential supersymmetric contribution, the lower bound on the mass of the lightest CP-even 
Higgs boson, and precision electroweak tests. We also require the lightest supersymmetric particle 
(LSP) to be the lightest neutralino. We do not, however, require that the thermal relic abundance tt x 
of the LSP falls within the CDM abundance determined within the ACDM paradigm. We assume that 
non-thermal production of neutralinos in the early Universe |35j , or cosmological enhancements of the 
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Figure 5: Sample soft (solid red lines, bb final state) and hard (dashed blue lines, t + t~ final state) photon 
spectra giving a photon flux compatible with the excess observed from the direction of Draco by the CACTUS ACT. 
The WIMP mass is set to 200 GeV, and the annihilation cross section times velocity to (o~v) = 10 -23 cm s . 
The dotted lines correspond to an estimate of the diffuse gamma-ray background, and the left panel refers to 
an angular resolution Afl — 10~ 3 sr, while the right panel to an angular resolution Aft = 10~ 5 sr, with the 
assumed NFW2 dark-matter halo. 





m\ 


m 2 


m 3 


rriA 


m s 




tan/3 


30 -r- 1200 


2 -T- 1200 


50 -r- 1200 


m LSP ~r~ 20000 


100 -T- 10m LSP 


(1 ~r~ 10)m LSP 


(—3 ~r~ 3)rrig 


1^-60 



Table 2: Ranges of the MSSM parameters used to generate the models shown in Figs. and0 All masses 
are in GeV, and m LSP = min(^, mi, 7712). The quantity indicates the following scalar masses (which were 
independently sampled): m Q 1 3 , m u 13 , m d 13 ' m L 12 3' m ei 2 3 ■ ^° av °id FCNC constraints, we assumed the 
squark soft supersymmetry breaking terms of the first two generations to be equal. Ag stands for the third 
generation sfermion trilinear terms: those of the first two generations were taken to vanish. 

relic neutralino density [HEj, brought f2 x ~ ^cdm- We detail in Table [2] the MSSM scan procedure. 

In the upper left panel of Fig- El we indicate on the (m x , {o~v)) plane supersymmetric models that 
give a gamma-ray flux from Draco in the range of Eq. (|12|) . Filled points are also consistent with the 
(strongest) requirement of a sufficiently low gamma-ray flux above 150 GeV. All models are consistent, 
instead, with the conservative constraint, from the Poisson variance of the measured counts, on the 
flux above 150 GeV. We employ here the BUR2 profile of Table ^ We also indicate the regions 
that significantly improve the fit to the gamma-ray flux from the direction of the Galactic center 



14 





OOOV 

oooooo&^ 
ooooooooodS 



OOOOOOOOOOOOOOOOOOO' 
ttOtO«OO0OOOOOOOOOOOO(?.. 

00«OOODOOOODDOOOi 

OOOOOOOOOOQC 

— otoooooooc 

O0OOOQ 



A 
> 
6 

V 



10 



EGRET preferred region (cored MW halo) 
EGRET preferred region (NFW MW halo) 



10" 



100 



Neutralino Mass (GeV) 



1000 



10" 



3 

o 

X io" ; 



a 

< 

> 

CD -! 

a io 

A 
W 



10" 





MOOOOOOOOOOOOO 


• •••••••••••••••oo»»oooooooo 

••••••••••••••••••oooooooo 

••••••••••••••••••oooooo 


GLAST sensitivity 


•••••••••••••••0O£ 

••••••••••••••c 

•o«*«*« 

OOOQ 

oo 





100 



Neutralino Mass (GeV) 



1000 



Figure 6: We show the results of a scan over the general MSSM (see Tabled for details) where we looked for 
supersymmetric models providing a flux of gamma rays from Draco compatible with the CACTUS excess. The 
left panel shows the (towimp, (o~v)) plane, while the right panel shows the expected integral flux of photons at 
GLAST, with the corresponding estimated point-source sensitivity in the direction of Draco. 

as measured by EGRET |37j . for a cored and an NFW Milky Way dark-matter halo. Intriguingly 
enough, for a Milky Way halo that is slightly cuspier than a Burkert (cored) profile (solid black line 
in Fig. H3 right), the region favored by the EGRET data can be consistent with the Draco signal, 
provided the neutralino mass is not much heavier than 250 GeV. We also indicate the upper limit on 
the annihilation cross section in the general MSSM derived in Ref. jS], and note that most models fall 
very close to the maximal cross sections in the context of supersymmetric dark matter. Such large 
cross sections could be in conflict with the production of other secondary annihilation products, such 
as antiprotons, positrons, or antideuterons |38U39| . However, the large uncertainties in the modeling of 
diffusion processes and nuclear reactions, together with those connected to the Milky Way dark-matter 
halo, can leave the freedom to circumvent those constraints |39| . 

We find that over the whole scanned supersymmetric parameter space, the total flux over the whole 
1 degree angular region features, with the BUR2 profile we employ here, an upper limit <j>% < 2x10- 10 , 
with a maximum at a neutralino mass of 600 GeV. This upper bound indicates that no supersymmetric 
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Figure 7: The neutralino -'proton scalar-interaction cross section and the flux of muons from the Sun induced 
by neutralino annihilations, with the corresponding present and future experimental sensitivities, for the same 
models as in Fig. 

models can give, in the present halo setup, 100% of the excess counts over the off-source background 
detected by CACTUS. The latter could be achieved with a cuspier NFW profile, the scaling among 
the fluxes simply being NFW /0 BUR = jNFW/jBUR ( gee pig.QJ. 

We then carry out an analysis of the prospects for the detection of the CACTUS-compatible 
supersymmetric models with other search avenues. In particular, in the right panel of Fig. we show 
the total photon flux expected at GLAST, integrated over an energy threshold of 1 GeV, and with 
Ail = 10~ 5 sr (which maximizes the signal to noise with the profile we use here). All CACTUS- 
compatible supersymmetric models will be, under the present assumptions, unambiguously within the 
sensitivity of GLAST. The upper bound on the integrated flux of photons stems from the EGRET 
bound, Eq. JT3J), since </> GLAST = 1O- 2 EGRET ( J_ 5 / J_ 3 ) ps 2.3 x 10~ 9 cnT^s -1 . 

Fig. [7| collects, again for the CACTUS-compatible supersymmetric models discussed above, the 
results for the rates at future direct-detection experiments (left) and at IceCube (right), where muons 
induced by neutrinos from dark-matter annihilation in the Sun will be sought. Although not guaran- 
teed, the detection prospects for supersymmetric models featuring the correct mass and annihilation 
cross sections certainly look promising. 

If we require the thermal neutralino relic abundance to agree with the CDM abundance, then 
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Figure 8: Scan of the general MSSM (see Table 0) showing the putative values of Jig needed to explain 
the CACTUS excess for super symmetric models featuring a thermal relic abundance fl x h 2 ~ f^cDM^ 2 - The 
horizontal lines indicate the values of J_5 for a sample of Draco halo profiles tltf - including the NFW2 and 
BUR1 profiles of Figs. andgj 

strong constraints are placed to the range of annihilation cross sections. One can then ask, for 
neutralino models with the correct thermal abundance, what values of J_5 are required to account 
for the CACTUS excess? We answer this question in Fig. El where we indicate, for MSSM models 
obtained in the same scan outlined in Table [21 that fulfill the condition £l x ~ ^cdMj the values of J_5 
such that <f)% = 2.4 x 10" 11 cm s . We show with vertical lines the values of J_5 for the profiles 
employed in Figs. |3 and 0J The vertical arrows indicate the overall ranges for J_5 obtained in the 
analysis of Sec. [2j Evidently, even the minimal gamma-ray flux compatible with the CACTUS excess 
requires an extremely cuspy halo for Draco. As we pointed out before, such cuspy profiles should be 
readily identifiable once a better angular resolution is achieved through the analysis of the angular 
distribution of the photon excess counts. 

3.4 Monochromatic gamma rays 

Given the preliminary spectral structure of the CACTUS Draco gamma-ray excess — i.e. a large amount 
of counts in a narrow energy range — it is certainly worthwhile to investigate whether the excess may 
be due to direct annihilation of dark- matter particles to photons. Unfortunately, the poor energy 
resolution of the solar-array ACT is not suitable for a prompt discrimination of this possibility, which, if 
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Figure 9: Left: The ratio between the number of photons from direct pair annihilation of neutralinos into two 
photons and from the gamma-ray continuum, integrated above 50 GeV, for the CACTUS- compatible supersym- 
metric models of Fig. Right: the ratio of the annihilation cross section of neutralinos into two photons over 
the total neutralino annihilation cross section, for the same ensemble of models. 

detected, would constitute a smoking gun for dark-matter annihilation. A dominantly monochromatic 
signal would moreover circumvent the friction with the EGRET null result. 

As a first step, we again consider here super symmetric dark matter, and determine whether a 
scenario where the bulk of the dark-matter-induced gamma-ray flux originates from direct neutralino 
annihilation into two photons is at all viable. To this extent, we show in the left panel of Fig. ED 
the relative fraction of monochromatic photons versus the integrated number of continuum photons 
with energies larger than 50 GeV, for super symmetric models giving a CACTUS-compatible gamma- 
ray flux, as determined in the previous Section. In the right panel, we show the relative branching 
fraction for neutralino annihilation into photons. We conclude that within supersymmetric models, 
the bulk of the photon flux always stems from the continuum component, the monochromatic part 
contributing less than ~ 0.2% of the photon counts for CACTUS-compatible supersymmetric models. 
As a byproduct, we derived from the full scan an upper bound on the branching ratio of neutralino 
annihilation into two photons in the general MSSM which reads 

r = (av) 77 /(av) t0t < 3 x 10" 4 . (16) 

On more general grounds, outside the supersymmetric paradigm, one can still hypothesize that the 
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Fi gure 10: Limits in the (jn x , (<7i>) 77 ) plane to the ratio r 



f /(crv} tot from the EGRET bound on the 



low-energy gamma-ray flux from Draco, for the cuspy NFW2 profile (left) and for the cored BUR1 profile (right), 
and for soft (solid lines) and hard (dashed lines) continuum photon spectra. 



bulk of the CACTUS signal comes from the monochromatic line, and constrain, through the EGRET 
bound, the continuum contribution, and hence the quantity r. We point out that the Z7 or the Zh 
lines would constitute less favorable scenarios here, since (1) the monochromatic photon flux would 
be smaller by a factor 0.5, and (2) the continuum photons from the Z and h decay would contribute 
to the continuum photon yield, strengthening the EGRET bound. 

We assume, in Fig. 1101 hard and soft photon spectra (dashed and solid lines: the regions lying 
above the lines are excluded by the EGRET constraint), and cuspy and cored profiles for the Draco 
dark-matter halo (left and right panels). The green regions correspond to the CACTUS signal: in 
the case of a cored profile, the branching ratio into two photons must be at least as large as 1-^10%, 
clearly incompatible with supersymmetric dark matter. Requiring that all of the CACTUS excess 
counts originate from photons produced in dark-matter annihilation, and consistency with the EGRET 
bound, implies much larger branching fractions (r ~ 0.5 — 0.99) and a dark-matter model where the 
dark-matter particle predominantly annihilates into two photons. 
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3.5 Decaying dark matter 

Uncertainties in the innermost structure of dark-matter halos critically impact the computation of 
the flux of photons from dark-matter annihilation. Since the annihilation rate per unit volume is 
proportional to the square of the dark-matter density, the occurrence of high-density spikes or cusps 
in the center of dark-matter halos may lead to flux enhancements of various orders of magnitude, 
without affecting significantly the outer structure of the dark halo, where the rotation curves are 
usually best determined. If, alternatively, the gamma-ray excess results from the decay of a quasi- 
stable dark-matter particle, the flux is only linearly proportional to the dark matter density (and 
hence, in the case of a distant object like Draco, to a first approximation, to the total dark-matter 
mass), which is much more reliably constrained by dynamical measurements, inducing a signficantly 
smaller spread in the flux predictions. 

We now show that the CACTUS Draco gamma-ray excess cannot be explained in terms of a 
decaying dark-matter particle. The resulting flux of photons in the diffuse gamma-ray background 
and from the center of our own Galaxy would in fact be exceedingly large. The ratio of the gamma-ray 
flux (jp raco from Draco's halo and that i^ MW from the Milky Way halo can be written as 

GC ^ B = W^ \w^)\L oi ^ p m ) ' (17) 

where the line-of-sight integral is performed either in the direction of the Galactic center (rcc) or 
in the direction of Draco, to evaluate the diffuse gamma-ray background (roB)- If the Draco mass 
is formed entirely by dark matter, then M^™ co ~ 8.6 x 10 7 M |25j . The distance to Draco is 
^Draco ^ ^ Q |4Q| ^he Mjiky Way line-of-sight integrals depend on the assumed dark-matter halo. 
We evaluated the integrals for two extreme choices, a NFW profile [2U] with scaling density and radius 
p s = 5.4 x 10 6 M /kpc 3 and r s = 21.8 kpc respectively, and a Burkert profile [2J2|, with scaling density 
and radius p s = 1.5 x 10 7 M Q /kpc 3 and r s = 11.7 kpc. We find 

r DB <1.5-1.7xlO- 4 , r GC < 0.25 - 68 x 1(T 6 , (18) 

with the smaller figures corresponding to the NFW Milky Way profile. Over a 1 degree angular region, 
the diffuse gamma-ray background from the direction of Draco quoted in Eq. © gives an integrated 
flux over 50 GeV 4>jyg < 10~ n cm~ 2 s -1 . We do not have a direct measurement of the gamma-ray 
spectrum from the Galactic center in the energy window between 10 GeV and 150 GeV, the lower 
limit being the highest energy probed by EGRET, and the upper limit the lowest energy probed by 
HESS. To have an idea of the integral flux above 50 GeV one can extrapolate the power-law behaviors 
featured by both the EGRET (at E > 2 GeV) and HESS differential photon flux, respectively, giving 

-3.1 



O*0egret -, a -in-6 ( E \ -2 —In \r—l 

~ 1.6 x 10 — — cm s GeV 



dE 



\GeVj 
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d<^HESS K in-12 ( E \ _2 _!„ , r _l 

~ 2.5 x 10 cm s lev 

d£ V T eV J 

These figures, extrapolated and integrated over the energy range E > 50 GeV, give the following 
integrated fluxes: 

Regret - 2 x 10" 10 cm-V 1 , 
0|ess - 8 x 10_11 cm^V 1 . 

The flux of gamma-ray photons from dark-matter decay from Draco is therefore conservatively bounded 
to be 

^Draco ~ 1.7 x cm^V 1 , diffuse gamma — ray background, (19) 

^Draco ^ 1A x 10_14 ™" 2 s _1 , Galactic center, (20) 

which is clearly incompatible with the estimates in Eq. ()12|) . 

We therefore conclude that the expected flux from the Galactic center and from the diffuse gamma- 
ray background are not consistent with the putative CACTUS signal from Draco in the context of a 
decaying dark-matter particle. A sufficiently large Draco flux would evidently violate the available 
constraints to the diffuse gamma-ray background and on the flux of gamma rays from the center of 
the Milky Way by various orders of magnitude. 



4 Conclusions 

We considered in this paper the possibility that the gamma-ray excess observed by CACTUS from 
the direction of the Draco dSph galaxy originates from WIMP annihilation. We summarize below the 
main results of the present analysis: 

• We showed that future measurements, with ~ 0.1 degree angular resolution, should allow us to 
distinguish between cored and cuspy halos, even though the current CACTUS angular resolution 
is not good enough. 

• We estimated the putative gamma-ray flux from dark-matter annihilation considering the two 
extreme possibilities that the dark-matter signal consists of (1) all the excess counts reported by 
CACTUS, and (2) the excess over background corresponding to Draco's innermost region. 

• We analyzed, in a model-independent approach, the regions of the particle-mass versus annihilation- 
cross-section parameter space compatible with the estimated CACTUS excess, imposing the con- 
straints from the null results of gamma-ray searches reported by EGRET and by the VERITAS 
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collaboration observation with the Whipple-lOm ACT, and from the absence of a statistically 
relevant excess in the CACTUS data for photon energies above 150 GeV. The annihilation cross 
section range lies well above that expected with standard thermal relic dark matter production. 

• The total excess counts over background reported by CACTUS can only be reproduced with a 
very cored dark-matter halo and a hard photon spectrum. If this is the correct explanation for 
the excess, then higher-angular-resolution measurements should still show a ~ 1 degree spread 
around Draco's center. Almost all CACTUS-compatible models will be within reach of GLAST. 

• In the case of super symmetric dark matter, the annihilation cross sections needed to reproduce 
the CACTUS signal are close to the maximal theoretically allowed values and would imply a 
negligible thermal relic abundance. CACTUS-compatible supersymmetric models give typically 
very large detection rates at direct-detection experiments, and a sizeable neutrino flux from 
neutralino annihilation in the Sun, which in principle could allow a cross-check of the gamma- 
ray signal that GLAST should see within this scenario. 

• The possibility of a dominantly monochromatic origin of the CACTUS excess is not viable 
within supersymmetry, and requires, for consistency with the EGRET null result from Draco, 
dark-matter models that annihilate predominantly to two photons and thus produce a very 
suppressed continuum photon spectrum. 

• A decaying-dark-matter scenario is ruled out by the resulting inferred gamma-ray flux from the 
center of the Milky Way and in the diffuse gamma-ray background. 
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